AD-A076  550  NAVAL  POST6RADUATE  SCHOOL  MONTEREY  CA  F/G  13/13 

THE  CONSISTENT  SECOND-ORDER  THEORY  OF  WAVE/STRUCTURE  1NTERACTIO— ETC (U) 
SEP  79  C  J  GARRISON  NSF-EN673-04019 

UNCLASSIFIED  NPS-69-79-010  NL 


1  OF  | 

AD 

AO 78660 

■ 

s'", 

AT 

* 

END 

DATE 

FILMED 

12  79 

DDC 

ADA076550 


©LEVEL 


•  NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THE  CONSISTENT  SECOND-ORDER  THEORY 


OF  WAVE/STRUCTURE  INTERACTION 


C.  J.  Garrison 


September  1979 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


Rear  Admiral  Tyler  F.  Dedman  J.  R.  Borstlng 

Superintendent  Provost 


THE  CONSISTENT  SECOND-ORDER  THEORY 
OF  WAVE/STRUCTURE  INTERACTION 

The  consistent  second-order  theory  of  the  Interaction  of 
regular  gravity  waves  with  a  fixed  object  in  water  of  finite 
depth  Is  developed.  The  theory  is  carried  out  for  the  most 
general  case  of  a  body  of  arbitrary  shape  which  may  extend 
through  the  free-surface  or  be  completely  Intnersed.  The 
Incident  wave  evolves  in  the  development  as  a  second-order 
Stokes'  wave.  Boundary-value  problems  are  established  for 
both  the  first-  and  second-order  velocity  potentials  and  a 
numerical  method  based  on  the  Green's  function  Is  outlined. 

The  work  reported  herein  has  been  oartiallv  supported  by  the 
National  Science  Foundation,  Washington,  D.C.  205S0. 


C .  Ga rr 1  son 

Associate  Professor  of 
Mechanical  Engineering 


Approved  by: 


JMar to.  Chairman 
Department  of 
Mechanical  Engineering 


Willi,*  W.  Tolies 
Dean  of  Research 


SCCOHITV  CL  1C  ATiOh  O*  This  p*ct  (»***  D«««  J; 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

»  «.!»:»•  n^mHah  i  govt  acc i show  wo 

IHO  NPS- 69  -79-31  P 

<  T|  T\_t  fmnO  TuA.Ift.) 

)  THE  CONSISTENT  SECOND-ORDER  THEORY  OF  WAVE/  1 

'  STRUCTURE  INTERACTION  , 

1  type  OF  HPOAT  4  PIKIOO  COVENO 

)  Final  rep«*V^Sept .  1 979 

1  *  , 

«.  a!a««awn«  oaa.  »e^w»r  wuwacn 

y  Au  Two*  •/ 

t  contract  on  caanT  nunICA.'O 

NSF  Grant  No. 
p-EN(C73-04019  ** 

t  AtWAOWWIWO  ONGAWIl  ATlOW  W«W(  ANO  AOOAftl 

Naval  Postgraduate  School 

Monterey,  CA  93940 

10  paogaam  KlCmint  PAO.^CT,  T Af k 
ARIA  i  *»0«K  UNIT  NUMBCMS 

/'  •  ' 

It  CONTROLLING  O'^lCf  AnO  A03«(St 

II 

National  Science  foundation  jj. 

Washington,  DC  20550 

U  4ET3AT  OATS 

i  ^ppMMMI7q 

a.  »acu  — - 

14  moniTONinG  AGInCt  N*wf  %  A30N  f  S  V  1  dtllf'ynt  twm  CoAfra/ffnj  Otttc*) 

14  ftCCJftlTV  CLASS,  cl  thl»  »•>•#<)  j 

Unclassified 

11  a.  0£CL  MSiHC  ATIOn/ODWnCNAOiNO 
icwi:<ki 

I*  OifT  m  (»..j  T|  on  UATtufuT  (•*  M  • 


Approved  for  public  release;  distribution  unlimited. 


17  *f  (*•  I-  !*/•<*  70.  1 1  *  #?«>'•, 


»•  Ho?n 


>t  *£v  »C*0i  *•••••  •  J»  •/  •*(*  I4»~fty  A*  4fac*  r  *+*••) 

wave/ body  interaction 
wave  forces 
nonlinear  wave  forces 
Green's  function 

10  AH|TN«Ct  I’Wtf*  (M  Il^t  If  ltf»  +  tlfy  if  blyrk  ngplar) 

The  consistent  second-order  theory  of  the  interaction  of  regular  gravity 
waves  with  a  fixed  object  in  water  of  finite  depth  is  developed.  The  theory 
is  carried  out  for  the  most  general  case  of  a  body  of  arbitrary  shape  which 
may  extend  through  the  free-surface  or  be  completely  immersed.  The  incident 
wave  evolves  in  the  development  as  a  second-order  Stokes'  wave.  Boundary-value 
problems  are  established  for  both  the  first-  and  second-order  velocity  poten¬ 
tials  and  a  numerical  method  based  on  the  Green's  function  is  outlined. 


1473  foi-os  o»  '  wow  o  ii  oasoitTt 

I / i  o  io  j-o  i  *•  »*fl  l 


Unclassified  -  ^  X-  L  . J 

»«CuAITY  CLAIliriCATICN  o*  Tni»  W ACC  D»*» 


DO  ,0”' 

v\j  |  JM  ?  | 


INTRODUCTION 

-"The  determination  of  forces  exerted  by  qravity  waves  on  large  structures 
immersed  in  the  sea  has  become  of  oreat  practical  interest  in  recent  years.  For 
example,  in  the  desiqn  of  bottom-mounted  oil  storage  facilities  or  larqe  ocean 
caissons,  the  wave-induced  horizontal  and  up-lift  forces  and  overturning  moments 
are  factors  of  primary  importance.  The  effect  of  larqe  amplitude  waves  in  parti¬ 
cular  is  of  importance  in  the  determination  of  the  permanence  of  an  ocean  struc¬ 


ture  and.  therefore,  a  hiqher-order  theory  appears  to  have  significant  practical 
value.  For  example.  Apelt  and  Mackniqht  (1976)  found  measured  forces  on  a  ocean 
caisson  model  in  fairly  large-ampl itude  shallow-water  waves  to  be  considerably 


in  excess  of  calculations  based  on  linear  diffraction  theory. 


With  this  application  in  mind  then,  the  solution  through  the  second-order 


in  wave  height  is  developed  herein  for  regular  wave  Interaction  with  a  fixed 


object  of  arbitrary  shape  in  water  of  finite  depth.  The  theoretical  development 
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is  carried  out  in  a  mathematically  consistent  manner  in  that  all  terms  through 
the  second-order  in  wave  height  are  included  in  the  analysis,  and,  as  such, 
represents  the  extension  of  the  linear  theory  presented  previously  by  Garrison 
and  Seetharama  Rao  (1971).  The  problem  is  recognized  as  a  regular  perturbation 
problem  in  the  small  parameter  (wave  height/characteristic  body  dimension)  and 
the  incident  wave  appears  within  this  framework  as  the  second-order  Stokes' 
wave.  A  numerical  method  for  solving  the  resulting  equations  suitable  for  digital 
computer  evaluation  is  outlined. 

REVIEW  OF  LITERATURE 

The  solution  of  the  linear  wave/structure  interaction  problem  is  now  fairly 
we  1 1 -developed  for  infinite  depth  as  well  as  for  the  finite  depth  case.  A 
number  of  papers  have  appeared  in  the  literature  dealing  with  the  oscillation  of 
two-dimensional  bodies  on  a  free  surface  in  water  of  infinite  depth,  a  problem 
which  is  ma thema t i ca 1 1 v  similar  to  the  fixed-body/wave  interaction  problem. 
Examples  include  the  work  of  Ursell  (1949),  Porter  (1960),  Vugts  (1968),  and 
Pauli  ing  and  Richardson  (1%2).  Dean  and  Ursell  (  1959)  made  both  an  experimental 
and  theoretical  study  of  small  amplitude  wave  interaction  with  a  fixed,  semi- 
submerged  circular  cylinder  in  deep  water.  Yu  and  Ursell  (1961)  treated  the 
problem  of  a  semi -submerged  circular  cylinder  oscillating  vertically  in  water  of 
finite  depth. 

Somewhat  less  attention  has  been  given  to  solvinq  the  corresponding  three- 
dimensional  problems.  Apparently  the  first  analysis  of  the  linear  wave  inter¬ 
action  with  a  fixed  body  was  carried  out  by  Havelock  1940)  for  the  case  of  a  pile 
in  water  of  infinite  depth  and  MacCamy  and  Fuchs  (1954)  solved  the  same  problem 
for  water  of  finite  depth.  Havelock  (1959)  also  evaluated  the  added  mass  and 
damping  coefficients  for  a  heavinq  semi  -  immersed  sphere  in  deep  water,  and  Wang 
(1966)  extended  this  to  include  the  finite  depth  case.  Haskind's  relations 
(Haskinds  (1957))  as  discussed  by  Newman  (1962)  may  be  used  to  determine  the 


linear  or  first-order  forces  actinq  on  the  fixed  body  from  a  knowledqe  of  the 
damping  coefficients  for  the  same  body  oscillating  in  still  water.  Garrison  (1974) 
has  made  this  evaluation  using  Havelock's  results  for  the  hemisphere  and  the  re¬ 
sults  were  in  good  agreement  with  his  direct  calculation  of  the  heave  force  based 
on  the  Green's  function  approach.  Kim  (1965,  1966)  has  applied  the  Green's 
function  approach  to  compute  the  excitation  forces,  added  mass  and  damping  co¬ 
efficients  for  a  semi -immersed  ellipsoid  in  deep  water. 

Most  work  on  wave/body  hydrodynamics  reported  in  the  literature  treats  con¬ 
figurations  and  conditions  which  have  primary  application  to  the  study  of  ship 
motion;  only  recently  an  interest  in  large,  fixed  offshore  structures  has  develop¬ 
ed.  Even  though,  probably  the  first  potential  flow  (or  diffraction)  solution  to 

1 

a  wave  force  problem  having  direct  application  to  large  ocean  structures  was  ob¬ 
tained  over  twenty  years  ago  by  MacCamy  and  Fuchs  (1954),  namely,  the  small 
amplitude  wave  interaction  with  a  fixed  vertical  circular  cylinder  (pile)  in 
water  of  finite  depth.  This  diffraction  solution  represents  the  only  closed  form 
solution  of  its  type  available  in  the  literature  and,  therefore,  is  of  primary 
importance  to  this  general  problem.  Little  work  of  this  type  having  application 
to  large  ocean  structures  has  appeared  since  that  time  until  recently  when 
Garrison  and  Seetharama  Rao  (1971)  applied  the  Green's  function  approach  to  calcul¬ 
ate  the  first-order  pressure  distribution  and  resulting  forces  acting  on  a  bottom- 
mounted  hemisphere.  Their  experimental  results  for  this  configuration  compared 
well  with  the  theory.  Garrison  and  Chow  (1972)  then  extended  this  analysis  to 
include  fixed  bodies  of  arbitrary  shape  and  compared  their  theoretical  results 
with  experimental  results  corresponding  to  two  different  submerged  oil  storage 
vessel  configurations.  Milgram  and  Halkyard  (1971)  also  have  developed  a  linear 
theory  and  applied  their  method  to  certain  axisymmetric  bodies  in  deep  water. 
Garrison  (1974,  1977)  has  more  recently  given  details  of  a  practical  method 
based  on  the  Green's  function  to  evaluate  added  mass  and  damping  coefficients  for 
floating  bodies  of  arbitrary  shape  in  water  of  finite  depth,  and  at  the  same  time 


calculate  the  linear  wave  forces. 


It  may  be  stated  that,  in  general,  linear  theory  is  quite  well  developed 
and  good  agreement  between  theory  and  experiment  has  been  obtained  for  wave  forces 
corresponding  to  small  amplitude  waves  in  both  the  two-  and  three-dimensional 
cases.  This  is  true  as  well  for  bodies  oscillating  in  a  free  surface  in  still 
water.  Experience  has  also  shown  that  linear  theory  gives  good  wave  force  results 
when  the  body  is  fairly  deeply  submerged  even  if  the  amplitude  of  the  incident 
wave  is  not  small.  However,  if  the  object  is  surface  piercing  or  not  deeply  sub¬ 
merged  and  the  water  is  shallow,  nonlinear  effects  caused  by  waves  of  finite  amp¬ 
litude  become  pronounced.  In  such  cases  linear  theory  is  inadequate,  and  a  higher 
order  wave/structure  interaction  theory  becomes  necessary. 

Nonlinear  potential  solutions  to  wave/structure  interaction  problems  are 
much  more  limited  than  linear  solutions.  Olgilvie  (1963)  solved  for  the  first- 
and  second-order  forces  on  a  horizontal,  submerqed  circular  cylinder  in  water  of 
infinite  depth,  the  wave  crests  being  parallel  to  the  cylinder  axis.  However, 
his  second-order  forces  included  only  the  time-average  part  and,  as  a  consequence, 
the  first-order  potential  only  was  needed  and  obtained.  This  work,  therefore,  did 
not  actually  represent  a  complete  second-order  solution.  However,  a  consistent 
second-order  theory  for  a  horizontal  cylinder  in  still  water  has  been  developed 
by  Lee(1968)  and  Potash( 1 970) .  Garrison  and  Smith(1977)  have  treated  the  corres¬ 
ponding  wave/fixed-body  interaction  problem. 

Recently  there  has  been  some  interest  in  developing  analyses  for  nonlinear 
wave  interaction  with  three  dimensional  fixed  bodies  in  water  of  finite  depth. 
Chakrabarti  (1972)  made  an  attempt  at  the  extension  of  MacCamy  and  Fuch's  (1964) 
linear  solution  for  the  vertical  circular  cylinder  (pile)  to  the  fifth  order. 
However,  he  treated  the  free-surface  boundary  condition  improperly  so  that  the 
results  are  mathematical ly  inconsistent  and,  consequently,  of  little  value  for 
making  judgements  regarding  the  magnitude  of  nonlinear  diffraction  effects.  More 
recently  Raman,  et.al.  (1976.  1977)  have  also  treated  this  problem 
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althouqh  the  method  of  solution  is  extremely  Comdex  and  there  appears  to  be  no 
simple  means  of  checking  the  results. 

Isaacson  (1977)  has  arqued  that  the  second-order  solution  for  a  body  with 
a  circular  waterline  curve  does  not  exist  because  the  free-surface  and  kinematic 
boundary  condition  on  the  cylinder  are  not  consistent  at  the  mean  waterline  on 
the  cylinder.  However,  this  inconsistency  does  not  preclude  the  existance  of  a 
solution;  it  appears  that  such  inconsistencies  at  the  point  where  two  different 
boundary  conditions  join  is  rather  common. 

In  this  report  the  solution  throuqh  the  second-order  in  wave  height  is 
developed  for  bodies  of  arbitrary  shape  held  fixed  in  regular  waves  in  water  of 
finite  depth.  The  problem  is  formulated  as  a  regular  perturbation  problem,  and 
the  incident  wave  is  shown  to  represent  a  second-order  Stoke’s  wave.  It  is  also 
shown  that  the  second-order  problem  is  the  same  type  of  boundary-val ue  problem 
as  the  first-order  problem,  the  primary  difference  being  that  the  second-order 
free  surface  boundary  condition  is  nonhomoqeneous .  The  solution  is  formulated 
for  both  the  first-  and  second-order  potentials  in  terms  of  a  Green's  function, 
and  the  numerical  results  are  presented  for  several  conf iourations . 


FORMULATION  OF  THE  PROBLEM 

The  problem  under  cons iderat ion  is  depicted  in  Fiqure  1.  A  riqid  object 
havinq  a  characteri stic  dimension  a  is  immersed  in  water  of  depth  h,  and 
a  train  of  regular  waves  propagates  in  the  positive  x-direction.  The  poten¬ 
tial  solution  is  sought  to  the  interaction  of  the  incident  wave  with 
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y  =  0  plane  mean  water  elevation 


r.URE  1  0EF IN  I T  IONS 


the  fixed,  rigid  body  through  the  second-order  in  the  wave  amplitude, 
or  more  precisely,  in  the  ratio  of  the  wave  amplitude  to  characteristic 
dimension  of  the  object,  S. 

The  fluid  is  assumed  to  be  irrotational  so  that  the  velocity 
potential,  «,  may  be  defined  as 

^  =  V<£ci,  g.  i  )  (j) 

where  q  denotes  the  fluid  velocity  vector.  Assuming  the  fluid  to  be 
incompressible,  it  follows  that  the  velocity  potential  must  satisfy 

V*4Xx(y  Z,  t  )  =  O  (2) 

within  the  fluid  region.  (The  barred  quantities  denote  dimensional 
quantities . ) 

On  the  bottom  the  kinematic  condition  specifying  zero  normal 
velocity  is 

•  */  *  >  ~  O  (3) 

where  h  denotes  the  mean  fluid  depth.  On  the  rigid,  wetted  surface  of 
the  body,  specified  by  S(i  y  z)  -  o  .  the  kinematic  boundary 

condition  describing  the  zero  normal  velocity  condition  is  given  by 

•  vS  =o  .  (4) 

The  elevation  of  the  free  surface  above  the  y  *  0  plane  is  denoted 
by  and  the  kinematic  boundary  condition. 


=  O 


(5) 


'  9  ]-  - 
Jy*l 

must  be  satisfied  on  y  ■  <i.  Moreover,  on  the  free  surface  the  pressure 
is  set  equal  to  zero,  and,  accordingly,  Bernoulli's  equation  supplies  the 
dynamic  free  surface  boundary  condition, 

<$.(*  l*,l)  +  i.l.if*  <$>,(?.  f, *,*)*+  (6) 


- * 

where  h  denotes  the  Bernoulli  constant. 

It  is  convenient  to  cast  the  boundary  value  problem  thus  far  established 
in  dimensionless  form,  and  for  this  purpose  Introduce  the  following  dimen¬ 
sionless  parameters: 

x  -  X/a  ,  y-  y  'a.  ,  *  s  'a  ,  <3  3/5  ,  h-h/5  - 

H-  h/zq.  ,  h*  ~  hT  a  ,  >'^aL3./q  ,  t  c rt  L  ^ 

<P  =  cr  ^  qd 


where  a  *  2»/T,  T  being  the  wave  period.  The  dimensionless  time  is  denoted 
by  t,  and  h*  denotes  a  dimensionless  Bernoulli  constant.  As  indicated  in 
Equation  (7),  all  the  coordinates  and  length  scales  are  made  dimensionless 
with  the  characteristic  dimension  of  the  fixed  object,  a. 

Using  these  parameters,  Eqs.  (2-6)  which  define  the  boundary  value 
problem  may  be  rewritten  concisely  in  dimensionless  form  as: 


.i,t )  - o 

in  the  fluid 

(8) 

<Pj*. >  •  o 

(9) 

* 


=  o 


on 


Scx,4lz;  =  o 


(10) 


where  e  relates  to  the  wave  amplitude  in  a  yet  undetermined  manner. 

It  may  be  noted  that  in  Eqs.  (11)  and  (12)  the  surface  elevation,  }  , 
and  consequently,  the  small  parameter,  t,  does  not  appear  explicitly. 
Accordingly,  it  is  necessary  to  expand  the  particular  functions  Involved 
so  as  to  form  functions  which  do  contain  c  explicitly.  For  example, 
substituting  Eq.  (14)  into  (13)  and  expanding  each  term  in  a  Taylor 
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series  about  c  ■  0  gives: 


•=■€  <$(*.o.l,i)  *~2l[  <%(*.£>.  I  t)  +  ?(*.  id)  <$  (*,3,1!,*  )]  t  Ou M 


(16) 


Similar  results  can  be  obtained  for  the  derivatives  of  *  which  appear  In 
Eqs.  (11 )  and  (12). 

Substituting  Eqs.  (13-16)  as  well  as  expressions  similar  to  Eq.  (16) 
for  the  derivatives  into  the  boundary-value  problem  given  In  Eqs.  (8-12), 
the  separate  boundary-val ue  problems  for  *2,  etc.  are  obtained.  The 
first-  and  second-order  problems  may  be  Isolated  and  are  given,  respectl vely, 


by: 

FI rst-order: 

vl4N*. y .  z  ,t )  -o  07) 

<p  ( X  -h  i  M-O  (18) 

'3  ' 

$(*  y  l  t)  ^  )  -  o  (19) 

t  )  -  >’  2  (  W  .  ■  =  o  (20) 

^  ;  o  (21) 

't 

Second-order : 

S7z4>(x  y  2  t  )  =r  O  (22) 

4>  ?  O  (23) 


JO 


=  O  on 


(24) 


(x.o.  e,n  -  \>  )ti*  i.t)  =  -  }(*  z.i )  <£<.  x.o.i.t )  +  <$(x  o.i,i)fi(x.i.{)  (25) 

'i  t  l  '»y  *  » 

*■  <p(*.o.  t.i )  ? 

+  4>(*  o,i  i)  =-)(*. 1 1  t.  <p(x,o,l,i)- (26) 

♦  <p(».o.  *,0*  /-  <$>(  r,o,l,t)Zl  v  >7* 

y  <  J 


Both  }(  and  are,  in  addition,  subject  to  a  suitable  radiation  condition 
which  limits  the  scattering  disturbance  to  outgoing  waves. 

A  great  deal  of  similarity  is  evident  between  the  first-  and  second- 
order  problems.  In  fact,  the  only  basic  difference  lies  in  the  right 
hand  side  of  the  free  surface  boundary  conditions;  in  contrast  to  the 
first-order  free  surface  boundary  conditions,  the  second-order  boundary 
conditions  are  not  homogeneous,  the  right  hand  sides  of  which  being 
dependent  on  the  first-order  potential. 

Before  proceeding  to  the  solution  of  the  boundary-val ue  problems 
developed,  advantage  may  be  taken  of  the  fact  that  the  first-order  poten¬ 
tial  will  be  represented  by  a  periodic  function.  Accordingly,  we  may 
define  the  complex  potential  u(x,y,z)  as 

<£(*  y,  i  .*)  -  a  b  Re[  i  u,  (x.y.-z)  0  *  ]  (27) 


where  Kp  denotes  the  real  part,  and  a  -  2«a/l,  a  being  the  characteristic 
dimension  of  the  body  and  l  the  wave  length.  The  symbol  b  denotes  an 
unknown  real  constant.  It  is  also  appropriate  and  comon  practice  In 
linear  interaction  problems  to  express  the  potential  as  the  sum. 


;  j 


cp  — 


(28) 


where  ♦*  denotes  the  incident  wave  potential,  and  ♦  *  denotes  the  scattering 
potential  which  is  due  to  the  presence  of  the  rigid  body.  The  complex 
potentials  u*  and  uS  are  then  defined  in  relation  to  i^and  **  according  to 
the  form  of  Eq.  (27)  so  that  it  follows  that 


u,  -  u 

1  / 


u 


(29) 


It  is  recognized,  moreover,  that  the  incident  wave  potential,  j1, ,  must 
satisfy  the  first-order  problem  when  no  body  is  present,  i.e.,  »x  must 
satisfy  Iqs.  (17),  (18),  (20),  and  (21).  These  equations  represent  simply 
the  boundary- val ue  problem  for  the  well-known,  first-order  progressive 
wave.  The  solution  to  tnis  in  terms  of  the  present  notation  is  given  by: 


Uf  =  -  —  e‘<)X  (30a) 

Q  Cosh (q  bj ) 

In  [q.  (30a)  a  is  defined  in  terns  of  v  by  the  familiar  expression  from 
1  Inear  wave  theory 


V  cr 1  3  /<0  =  a  fanhcah) 


(30b) 


where  h  denotes  the  dimensionless  water  depth. 

Now,  substituting  Eqs.  (30a),  (29),  and  (27)  into  the  first-order 
problem  defined  by  Eqs.  (17-21),  and  eliminating  n  between  Eqs.  (20)  and 
(21),  the  following  boundary-value  problem  is  established  for  the  first- 
order  complex  scattering  potential: 


1 2 


V1u'c*.4,i  )  =  o 


(31) 


obtain,  upon  eliminating  n  between  tqs.  (25)  and  (26),  the  following 
boundary-value  problem  for  the  second-order  incident  wave  potential: 


z.t  >  -  O  (36) 

4?  (x,-h  z  t  )  o  (37) 

v 

<4>\*.o.  /  t)  t  v*  <p  («.o,  t  ,t)  =  -  4  (a1  -  vl)  smlziax-  t  >]  (38) 

*1  ltt  * 

The  periodic  solution  to  the  boundary-value  problem  specified  by 
tqs.  (36-38)  may  be  expressed  as 

<t>  =  j.  a  &  y>  r?e/  i  e~Ui  j  (39) 

where  the  second-order  complex  potential,  u*, ,  is  given  by 

a**  <i)  =  -  -i-.  cpsMzac'i+v)}  (40) 

smhcah) 

This  expression  is  familiar  in  wave  theory  and  is  of  exactly  the  same 
form  as  that  given  by  Lagleson  and  Dean  (1966)  for  second-order  Stokes’ 
waves . 

To  complete  the  establishment  of  the  boundary-value  problem  for  the 
second-order  scattering  potential,  it  is  necessary  to  substitute  the 
forms  j  *  i  *  and  t ,  *  ♦  «s  into  tqs.  (22-26)  and  eliminate  n 

l  *  •  *■  2  2  2 

between  tqs.  (25)  and  (26).  In  addition,  the  known  expressions  for  the 
incident  wave  potentials,  tqs.  (30)  with  Eqs.  (29)  and  (27)  for  the  first- 
order  potentials,  and  for  the  second-order  potentials,  tqs.  (40)  with 
tqs.  (39)  and  (35),  are  utilized  to  obtain: 


Vl<P*(x.  H  l-t  )  =  O 


(41) 


<p  (x  -h  l ,  t  )  •  O 


(42) 


<£*(», u  z  t)=--a.i2i^Jl^lj'^.coshI/a^MJ^L^,s.nhLztt(MjeUl*,  ^1(43) 

'  *  S»nh*con)  “• 


<# 


r  y<$  it  <4*  ^  *-  «Y  "  5  v  ^'y  ) 


vt  <e<  <  <  -*< 


(44) 


U*  -  3  u. 
i  y 


!* )  eui]]  +  U(k.b) 


In  £q.  (44),  U(x,z)  Is  a  time  independent  function  which  is  generated  by 
substitution  of  t ;  into  Eqs.  (2b)  and  (26),  However,  since  the  function 
is  not  needed  in  the  second-order  theory  for  purposes  of  evaluating  the 
pressures  and  forces,  It  is  for  brevity’s  sake  not  written  out  In  full 
here. 

The  boundary-val ue  problem  now  established  for  as  given  in  Eqs. 
(41-44)  is  linear,  and  except  for  U(x,z)  occurring  In  Eq.  (44),  is  time 
dependent  like  e'12t.  Accordingly,  It  is  appropriate  to  express  the 
solution  in  the  form 

<p  r  d  Q  6‘  y  Rcji  *■  i  (45) 

where  the  second  term  denotes  the  part  of  the  complex  potential  which 
Is  Independent  of  time.  By  substitution  of  Eq.  (45)  Into  the  boundary- 
value  problem  described  by  Eq.  (41-44),  separate  boundary-val ue  problems 
arise  for  \i]  and  u*  .  However,  as  will  be  evident  subsequently,  ♦* 


only  is  needed  in  order  to  evaluate  the  hydrodynamic  pressures  and 
resulting  forces  to  the  second  order  and,  consequently,  the  tine 
independent  part  of  **  is  of  no  interest.  We,  therefore,  dismiss  it 
from  further  consideration  and  concentrate  on  the  solution  us. 

Substituting  Lq.  (45)  into  Fq.  (41-44)  the  follov/ing  boundary 
value  problem  is  established  for  the  second-order  scattering  potential: 


y‘u'(vl47)  -.o  (40 ) 

U*  (x,-h, *  >  =  O  (47  ) 

ul  U.H.i)  =  — Xx - j  H  SinhUa(K*y)]  *  i  n,  ooshUuCK,  y)  M  Za  *  (4P,) 

n  Siangan)  1  J 

I  # 

Uj^K'O, /)  -  4  V  u/(*  o,i)  r  ( x.o  (49) 


where 


4  a[;,  U.;  U‘n  .  u,'(u:v  -svu.\).  -L  «;  (  u.‘n-»uf,)  (so) 

-  .?  a!)  -  e  a?)  -  3  u,’*  J,.  „ 


-  4  6/,  U 


The  radiation  condition  limits  the  scattered  waves  to  outgoing 
regular  waves  and  is  given  by 


where  a^,  denotes  the  dimensionless  wave  number  for  the  second-order 
probliv  which  is  uefined  by 


(w) 


INCIDENT  WAVE 


At  this  juncture  it  is  appropriate  to  completely  define  the  incident 
wave  and  specify  the  unknown  constants  b  and  h*.  Solving  Eqs.  (21)  and 
(26)  for  nj  and  n2,  respectively,  and  substituting  the  results  into  the 
expression  for  the  elevation  of  the  free  surface  defined  by  Eq.  (M) 
gives 

O  -  "  £  1,1)  *  C‘  J-  4>C*.Oliti)  * 

(53) 

-J-  t  <Px>.o.i ,{ <$<■  r.°-  l.t)  J+ 


Equation  (S')  expresses  the  free  surface  elevation  in  terms  of  the  total 
potentials  and,  consequently,  includes  the  effect  of  the  scattered  wave 
as  well  as  the  incident  wave.  However,  it  is  presently  of  interest  to 
obtain  an  expression  for  the  free  surface  elevation  of  the  incident  wave 


alone  in  the  absence  of  the  object,  for  this  purpose,  therefore,  we  set 

$S  =  =  0  and  evaluate  Eq.  (53;  using  the  known  expressions  for  the 

■  * 

incident  wave  potentials,  Eqs.  (30)  with  (.'7)  and  Eq.  (40)  with  (39). 


The  result  is: 


,  Cfthlah)  U^zhltahl  Rcrc,’<‘>-  l,j }  .  0«  • 
4  &ir\H  (ah)  J 


where  n1  denotes  the  dimensionless  surface  elevation  due  to  the  inci¬ 
dent  wave  in  the  absence  of  the  body.  However,  the  constant  term  In 
Eq.  (54 )  must  vanish  so 


to* 


b*  c  a1,  -  y1; 

4-y 


(55) 


The  remaining  terms  in  Eq.  (54)  represent  sinusoidal  variations,  the 
second-order  contribution  navinq  twice  the  frequency  of  the  first. 

If  we  define  the  wave  height  as  the  elevation  difference  between 
the  trough  and  crest  of  the  incident  wave,  then  cb  must  represent  the 
wave  amplitude  since  the  second-order  contribution  to  the  surface 
elevation  in  Eq.  (54  )  is  the  same  at  both  the  crest  and  trough,  the 
difference  being  zero.  Therefore,  in  terms  of  the  dimensionless  wave 
height  as  defined  in  Eq.  (7)  we  have: 

£  b  =  H  =  H  / ZOl  (55) 


wnere  H  denotes  the  elevation  difference  between  the  crest  and  trouqh. 
Tne  incident  wave  profile  is  then  specified  by  the  dimensionless 
surface  elevation  as: 

t)  *  m  cos(a»-  t )  ♦  Ml  -A  Cx>^t[{ah)  +  co sh(zah))  (57) 

4  sm h\ah)  v 

.  cosUca*-i>j  +  OCh*) 

Equation  (57)  agrees  with  the  expression  for  the  second-order  Stokes' 
wave  given  by,  for  example,  Eagleson  and  Oean  (1966). 


PRESSURE ,  FURCES  AND  MOMENT 


The  forces  and  moments  acting  on  the  Immersed  surface  are  deter¬ 
mined  by  carrying  out  surface  integrals  of  the  pressure.  For  this 
purpose  the  pressure  is  obtained  from  Bernoulli's  equation  as 

=  -/><£-  <&*  *■  *  4>i  ]~  pgy  *  pgh  *  (58) 

This  expression  may  be  cast  in  dimensionless  form  by  use  of  Eq.  (7), 
(13),  (14),  (15),  (27),  (35),  (45),  and  £6),  and  the  result  carried 
to  the  second-order  in  wave  height.  The  resulting  expression  is  then 


>(*,  u.l  t  )  -  - 


y  -  Pep.  <  -< 


z 

l 


*LX 

a 1 


-  i 
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The  dimensionless  pressure  coefficient  expressed  by  Eq.  (59)  is 
defined  as 


r 


p 

pga 


(60) 


in  which  P  denotes  the  fluid  pressure. 

The  first  tern  in  Lq.  (59)  denotes  the  hydrostatic  pressure  as 
y  represents  tie  dimensionless  depth  beneath  the  mean  free  surface. 

Tne  second  and  third  terms  are  harmonic,  the  second-order  part  having 
twice  the  frequency  of  the  first-order  part.  The  final  second-order 
term  in  Eq.  (5Q)  is  independent  of  time  and  gives  rise  to  time-average 
or  steady  state  components  of  force  and  moments. 

The  components  of  the  dimensionless  force  and  the  moment  vectors 
may  now  be  expressed  in  terms  of  integrals  of  the  pressure  over  the 
wetted  surface  area.  That  is,  we  may  write 


(59) 
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c/t)  =  ff  r  3.  ds  ,  i  •  *.*. 


where  the  dimensionless  force  coefficients  associated  with  the 
x,  y,  and  z  directions  are  defined,  respectively,  as 


PS  <5 


(62a) 


PS 


(62b) 


-JkJLL 
P  S  a* 


(63c) 


in  whicn  F^(t),  Fy ( t ) ,  and  F  z ( t )  denote  the  tnree  components  of  force. 
Similarly,  the  moment  coefficients  are  defined  as 


c  a)  - 
L*'1 


(62d ) 


M*(t) 

pgQ4 


( 62  e ) 


Cji)- 


A 1t(t) 

pqs* 


(62  f) 


in  which  Mx(t),  My(t),  and  Mz(t)  denote  the  three  components  of  the 


wave-induced  moment. 
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The  functions  are  defined  as 


=  n«  ,  ‘  "*  (**) 

^  .=  Cd  ♦  y )  n,  -  £  ,  g^zn.-^n,  ,  xn,-(^)n,  (63b) 

The  forces  and  moments  In  Eq.  (6?)  are  defined  relative  to  a  coordinate 

system  parallel  to  0(x,y,z)  but  shifted  downward  a  distance  d  as  Is 

evident  from  Lq.  (63b).  The  unit  vector  directed  normal  to  the  wetted 

♦  ^ 

surface  is  defined  by  ft  -  i  *7,  ♦  ■'  ^  *  *  na  and  dS  denotes  an 
element  of  surface  area  made  dimensionless  with  the  characteristic 
dimension  of  the  body  squared,  a7. 

We  may  make  an  initial  ordering  of  the  force  or  moruent  coefficient 
as  given  in  tq.  (bl)  in  terms  of  tne  small  parameter,  c,  or  equivalent 
H,  as 

c."’--//5  r.  % JS  -  Hf)sr  %ds -*'/( n  rds-oo<‘>  (M) 

where  the  pressure  coefficient  is  expressed  as 

-r>  --  *  h  f-'  *  Hl  /;  *  H  (H*)  (6b ) 

/  i°  / '  / 

and  the  coefficients,  p0,  p(,  and  p?,  are  defined  in  an  obvious  manner 
by  comparison  of  Lq.  (6b)  with  Eq.  (39). 

It  may  be  noted  at  this  point  that,  if  the  body  is  completely 
submerged,  the  surface  integrals  in  Eq.  (64  ;  refer  simply  to  the  actual 
surface  area  of  the  body  which  is,  of  course,  specified.  However,  If 

21 


the  object  is  surface  piercing,  S  denotes  the  actual  wetted  area  which 
is  dependent  on  the  instantaneous  surface  elevation  of  the  water  on  the 
body.  Accordingly,  it  is  necessary  to  express  the  limits  of  Integration 

on  the  wetted  surface  S  in  terms  of  n  and  carry  out  the  expansion  retain- 

•> 

ing  terms  with  coefficients  up  through  H  . 

For  this  purpose  the  dimensionless  differential  surface  area  may  be 
written  as  dS  ■  dc  dl  where  dc  denotes  a 

dinensi onless  differential  arc  length  on  the  wetted  surface  in  the 
horizontal  plane,  and  dl  denotes  a  dimensionless  differential  length 
in  an  orthogonal  direction.  Thus,  it  is  possible  to 

express  dl  in  terms  of  y  so  that 

d  S  dc  df  =  dy  (56  ) 

V  /  - 

The  integrals  indicated  in  Lq.  (M  )  are  now  carried  out  over  the  wetted 
surface  with  y  running  from  y  *  -e  to  y  *  n.  Accordingly,  a  typical 
integral  in  tq.  (h*)  may  be  expressed  in  the  form 


*  -  -  e 


where  e  denotes  the  lower  limit  of  the  immersed  surface. 

The  integrals  of  the  type  given  in  Lq.  (67)  are  next  expanded  in  a 
Taylor  series  about  M  ■  0  using  Leibniz's  formula  for  purposes  of 
evaluating  the  derivatives  of  the  integrals.  Carrying  out  these 

expansions,  substituting  the  results  into  Eq.  (65  ),  and  retaining  terms 

2 

up  to  order  H  yields: 
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where  S0  denotes  the  surface  area  below  the  plane  y  ■  0  and  C0  denotes 
the  closed  waterline  curve  formed  by  the  intersection  of  the  y  ■  0 
plane  with  the  surface  of  the  rigid  object. 

The  first  term  in  Eq.  (•  )  represents  the  hydrostatic  force  or 

moment  on  tne  body  where  the  associated  displaced  volume  Is  defined  as 
that  beneath  the  y  «  0  plane.  However,  since  our  Interest  Is  in  the 
dynamic  forces  and  moments,  we  may  disregard  this  first  term  with  the 
understanding  that  tne  buoyant  force  defined  as  such  should  be  accounted 
for  in  order  to  determine  the  final  total  force  or  moment. 

Finally,  applying  Eq.  (21,  and  (27)  as  well  as  the  definitions 
for  p0,  pt  and  p2  we  obtain: 


L  (t)  = 


-  O,  ’  6 1.\  )<T  '**  !  Q  dS  *■  a1  f\ L _ ,  / 

T  1  *  2  Jc[e}rrny'  'J 

.  R.[u‘C-“]9:  dc  [,-y^r  -  i]  VM*  9.dc 

*  4^7  jj  \a't\L*  g  -  1J9, ds) 


Vs) 


+  O(h') 


2  J 


Equation  (69a)  may  also  be  written  in  coefficient  form  as 


CM)  = 


C£>sCS,t  l  )  f  H'fc  CVS (6]  -  +  ff* y  r-  Qt**) 


(69  b) 


where  the  first-  and  second-order  force  (or  moment)  coefficients  and 
phase  shift  angles  are  defined  by  comparison  of  Eg.  £9  a)  with  b) 
as  follows: 


(69c) 


wnere  the  dimensionless  force  coefficients,  Fj  and  F2^,  are  real. 

The  steady  state  (nonperiodic)  force  (or  moment)  coefficient  is  also  defined 
by  comparison  of  Eqs.  ('  ‘a)  and  (69b)  as 


fT" 


llGvhr  -  0  lu‘i‘ dc 


(69  e) 


The  first  term  in  Lq.  £9b)  represents  the  linear  solution  having 
tne  fundamental  frequency,  o.  The  next  term  which  represents  a  second- 
order  contribution  to  the  force  or  moment  Is  also  harmonic  but  at  twice 
the  fundamental  frequency.  The  last  second-order  term  represents  the 
steady  state  contribution  which  is  independent  of  time.  This  latter 


force  is  generally  referred  to  as  a  drift  force  in  ship  hydrodynamics. 
(See,  for  example,  Maruo  (I960).)  However,  drift  force  is  generally 
calculated  by  applying  the  momentum  equation  to  the  diffracted  wave 
rather  than  by  using  near  field  results. 

The  evaluation  of  the  force  or  moment  coefficients  defined  by  Eq. 
(5?c-e)  represents  the  primary  object  of  this  paper.  However,  in 
order  to  carry  out  the  evaluations  indicated,  it  is  clear  that,  given 
the  geometry  of  the  rigid  immersed  object,  it  is  necessary  to  evaluate 
u4  and  u2  as  well  as  the  derivatives  of  uj  on  the  surface  area  denoted 
by  SQ  and  along  the  waterline  curve  CQ. 

SOLUTION  TO  THE  FIRST -ORDER  PROBLEM 

Having  now  established  the  need  for  the  potentials,  u,  and  u;,  we 
return  to  the  consideration  of  the  solution  to  the  boundary-value  prob¬ 
lem  developed.  The  practical  method  of  solution  of  the  first-order 
problem  given  in  Eqs.  (31-34)  is  discussed  in  depth  by  Garrison(1974,1978) 
and,  therefore,  only  the  major  steps  will  be  outlined  here.  Following 
the  Green's  function  method  of  solution  we  write  u*,  as  the  Integral 
over  the  i-nersed  surface  SQ  as: 

ij.l)  -  JL.  f(  ^  (70) 

4fr  JJ% 

where  U.n.c)  denotes  points  on  the  irmersed  surface,  fx  (c.n.t)  denotes 
the  unknown  source  strength  and  dS  *  dS/a*  denotes  the  differential 
surface  area  on  the  ironersed  surface  made  dimensionless  with  the  charac¬ 
teristic  dimension  of  the  object,  a.  The  function,  G  (x,y,j;C,n,t'.v), 
denotes  the  Green's  function  which  must  satisfy  the  equation 


^  G<*. y. *;!,!.  s;  >0  Sc»-§>  6<«- 1)  (7!) 

in  which  5  denotes  the  Oerac  delta  function.  The  function,  G(x,y,*;t,n,;*,v) 
nust  also  satisfy  £qs.  (32),  (34),  and  the  radiation  condition.  Such  a 
function  Is  given  by  Wehausen  and  Laltone  (1960)  as 


j;  v>  =  -y  *■ 


*  v] ‘ Z(~*>  d«  (72.) 

J  -  v  ccsnc«h) 

+  l  zrr  ‘  Q'~  )  c°shla-  co  s*[Q.(**hi]  j  (an) 

Q'h  -  Vl/7  +  y> 


in  which 

R  =[(*-  */  *  C  y-//  *  (i  -  S)1] ,/L 


(7?b) 


R  =  +  (Z-Z  )lj 
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(7?c) 


*  =  [(*-  i/  *  (a  If]"1 


(72d) 


The  symbol,  alt  Is  defined  In  terms  of  h  and  v  as  the  solution 


2* 


I 


of  the  equation 

a,  ianh (Q^h )  -  V  e  o  (72c) 

In  view  of  Eqs.  (7  e,'  and  (30b)  the  symbol  a4  is  clearly  equivalent 
to  a.  However,  in  the  cast  of  the  second-order  Green's  function,  the 
equivalence  docs  not  hold  and,  therefore,  separate  notation  is  maintained. 
In  Eq.  (7J  a),  P.  V.  denotes  principal  value  of  the  integral. 

An  alternate  series  form  of  the  Green's  function  is  also  given  by 
Wehausen  and  Laitone  (1960)  as 

G('  7  J;  v;=  ZZ&'l-Qt'j&hliLi  £•***.# /».*;  -  j  Jm(a.n)l 

'  nxh  _  c>‘h  *  v>  L  °  J 

CK> 

*  +  Z 

A  »/ 

where  J0  and  denote,  respectively,  Bessel  functions  of  the  first  and 
second  kind  of  order  zero,  and  KQ  denotes  the  modified  Bessel  function 
of  the  second  kind  of  order  zero.  In  Eq.  (73a)  aj  is  defined  by  Eq. 

(7?e)  and  the  cuanti ties.x^  are  defined  as  real  positive  roots  of  the 
equation 


far?  (  M*  h)  f-  V>  =■  Q 


I 

(73b) 


The  solution  to  the  boundary-value  problem  stated  in  fqs. (31-34) 
is  given  by  fq.  (/"*);  it  remains,  however,  to  determine  the  source  strength 


il 


I 


function  in  order  to  evaluate  the  potential.  This  Is  accomplished  by 
application  of  the  boundary  condition,  Lq.  (33),  which  results  in  the 
following  integral  equation  from  which  f,  may  be  determined: 


J-ff  (if,  nLH' 


)  «  »/  «'*/'/  T;y)dS  ■  smhLachty)]* * n.  cashra^l  ( , 

•  Ccsn(ah)L  •> 


The  normal  derivative  of  0  evaluated  on  the  imersed  surface,  as  required 
in  tq.  (7*),  is  determined  in  a  strai ghtforvard  manner  by  differentiation 
of  either  Tq.  (7a)  or  (7a). 


solut:  •.  to  si  j o:*r  ordlr  pROt  Lcn 

The  solution  to  the  second-order  problem  may  he  formulated  by  appl ication 
of  Green's  theorem  to  the  fluid  region  denoted  by  Y  in  figure  2.  Applying 


FIGURf  2  Region  for  application  of  Green's  theorem. 

Green's  theorem  with  subjects  u^,  and  G,  (the  second-order  Green's  function) 
and  assuming  that  G,  satisfies  the  equation, 

V’Cit(x,y,  Z;5,1,S)  -  &[x-s)  (75) 

as  well  as  the  kine  it  ir  condition  on  V  given  by  Lq.(47),  the  radiation 


Pi) 
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condition  qiven  by  Eq.(51)  on  S^.  and  the  homogeneous  free  surface  con¬ 
dition, 

6m/*,*,*; 5,T,S>  -  O  (76) 


gives  the  result: 


J-[  l ufojtT)  UT)-  <£ «V, *, 5. 1 7)djU  J&S 

^  (77) 

-  -l J~{4{f,c,r)  -  *//s, *,j)]GJ*.v, i)  c/s 


The  first  term  in  Eq . (77)  may  be  considered  to  represent  a  source  and 
doublet  distribution  over  the  immersed  surface  of  the  body.  However, 

Lamb  (1932)  has  shown  for  the  case  of  an  infinite  fluid  that  the  potential 
mav  be  represented  as  either  a  source  distribution  only,  a  doublet  distri¬ 
bution  only,  or  sow  combination  of  both.  This  proof  is  easily  extended  to 
the  present  case  allowing  the  first  term  to  be  written  as  a  distribution  of 
sources  only.  Then,  if  the  free  surface  boundary  condition,  Eq.(49),  is 
used  in  the  second  integral  the  result  is 

tf«Vr.y,0  -  /r  j  {JUT)  <2<(*.V ds 

(78) 

-+J  f  ‘ft,  f)  ds 

*3 

in  which  is  the  unknown  source  strenath  distribution. 

The  second-order  scatter  ootential  defined  bv  Eq . ( 78 )  satisfies  Eqs. 

(46,  47.  49-51).  The  remaining  kinematic  boundary  condition  on  the  surface 
of  the  body  results  in  the  integral  equation 
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^  ff  I fj;  Wi) 


Sinh4(Qh) 


[fly  fl*  CaslfttfJ'+'fijJ*- 


(79) 


H4X 


which  is  to  be  satisfied  for  x.y.z  on  S,  . 

In  both  Eg . ( 73)  and  (79)  a  practical  interpretation  of  the  limits  of 
the  free  surface  area,  £s  .  should  be  understood.  In  deriving  Eq.(78) 
through  application  of  the  Green's  theorem  the  outer  limit  of  S$  was 
considered  to  be  larqe  enouqh  that  the  radiation  condition  was  satisfied 
on  .  In  practical  application  this  occurs  rather  rapidly. 

In  view  of  the  conditions  which  were  placed  on  the  Green's  function  in 
the  above  development  G-  is  clearly  verv  similar  to  the  Green's  function 

C 

used  in  formulatinq  the  first-order  problem.  In  fact,  if  V  is  replaced 
by  4v*  and  is  replaced  bv  QL  in  either  E q .  ( 72 )  or  (73)  the  re¬ 

sulting  function  is  the  appropriate  Green's  function  for  solution  of  the 
second-order  problem.  That  is. 

(so) 


where  G  is  qiven  by  Eg. (72)  and  (73)  and  <3tis  defined  through  Eg.(7?e) 

at  Jun/t  QLh  -  4  V  *  O  (81a) 

In  the  case  of  the  alternate  series  form  of  given  in  Eq.(73),  the  roots 
are  defined  throuqh  Eg. (73b)  as 

*Mk  Jun  h )  V-  —  O  (81b) 
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An  alternate  form  of  Fq.(7?a)  for  the  sources  on  free  surfaces  Is  given  by 

3,4V,  a.)*  z  ryj X<g+L 4* 

(82) 

.  t  27T  CU  CoshfatCytk)]  smk(**k)  JmHua) 

4-v/r  +■  snr\hx’axh 

This  if  the  form  of  the  Cireen's  function  which  is  obtained  by  setting  J  =■  c 
in  Eq.(72a)  and  v/hich  was  also  given  by  Wehausen  and  Laitone  (1960)  for 
the  case  of  a  harmonic  pressure  distribution  on  the  free  surface.  The 
nonhomoqeneous  form  of  the  free  surface  boundary  condition  given  in  Eq.(49) 
is  equivalent  to  the  boundary  condition  for  the  first-order  problem  for  a 
pressure  distribution  of  the  free  surface. 


numerical  PROCEDURE 

Eq.(69)  clearly  indicates  that  the  determination  of  the  forces  and 
moment  acting  on  the  immersed  obgect  rests  on  the  determination  of  the 
potentials  Up  i,  and  derivatives  of  u^  at  points  on  the  immersed  surface. 
The  first-order  scatterino  potential  is  specified  by  Eq.(70)  in  terms  of  the 
first-order  source  strength  function,  fp  Thus,  to  determine  u|,  it  is 
first  necessary  to  solve  the  integral  equation,  Eq.(74),  for  f.. 

We  may  develon  a  numerical  solution  to  Eg. (74)  beginning  with  the 
partitioning  of  the  immersed  surface  into  N  subdivisions  or  facets  of 
area  ASj  .  each  with  a  nodal  point  at  its  center  located  at  the  point 
(xp  y.,  /.),  Recognizing,  moreover,  that  f^  (f,J,J)  is  a  well-behaved 
function  for  smooth  bodies  we  may  define 


(83) 


hi  -  _L_  /V?y  <*t  ,  y, ,  l.)  l  A(t,  «*)  Ot>shfr(l*y^ 


(LX< 

(84) 


» 


fi.  =  fie yj ,  % ,  ij ) 


and  accordingly  we  ray  approximate  E q . (74)  by  the  complex  matrix  equation: 

(XLj(  y)  ftj  -  (8r,) 

Once  has  been  evaluated  by  use  of  Eq.(83)  the  solution  of  Eq.(8G) 

may  be  carried  out  to  determine  fj  •  at  each  nodal  point  on  the  surface  of 
the  body.  A  discussion  of  the  details  of  the  numerical  evaluation  of  /"\ 
including  the  evaluation  of  the  i/R-  singularity,  has  been  given  by 

Garrison  (197?.). 

Having  determined  the  source  strength  function,  the  first-order  poten¬ 
tial  and  its  derivatives  on  the  surface  of  the  object  may  be  determined  from 
Eq.(70).  For  this  purpose  we  may  replace  the  surface  integrals  with 


summations,  writing 

a,s  =  p  o')  £  . 

4  •  ».  j 


1 ,  j  - '  t. 


<.  =  f.  -  '• 


j  *•  >.  /,  /v 


uf  =  &  fy)  f  ,  i 
'  y<,  '■> 


■>  -•  >.  f,  /v 


U‘it  c  fi 


t  J  *  >,  t,  »  •  •  A/ 


where  u^  ,  u*  ,  etc.  denote  functions  evaluated  at  the  i*^  nodal  point 

1  K, 

on  the  immersed  surface,  S0.  The  complex  matrices  occurring  in  Eqs. 

Cl7  -90  )  arc  given  by 


/Jas 

J 

y  $pff  'is 


I 


/Vy; = 


,  i.  ;u,  T;  y)  dS 


(93) 


(9/|) 


In  evaluation  of  the  integrals  in  E qs . ( 83 )  and  (91-98)  it  is  generally 
adequate  to  simply  evaluate  G  and  its  derivatives  at  the  centroid  of  the 
panel  and  multiply  by  its  area.  However,  when  point  i  is  either  equal  to 
j  or  near  j  it  is  necessary  to  take  more  care  in  integrating  the  1/R- 
singularity  in  G.  The  procedures  given  by  Garrison  (1978)  were  used  to 
evaluate  the  integrals  on  quadralateral  panels. 

To  obtain  a  solution  to  the  second-order  problem  it  is  necessary  to 

* 

first  evaluate  the  function  f  (J,3)  required  in  Eq.(49)  and  as  defined  by 
fq.(50).  For  this  purpose  we  proceed  again  numerically  by  dividing  the 
mean  free  surface  (or  y  =  0  plane)  in  the  vacinity  of  the  body  into  K  area 
elements  and  evaluate  f  at  the  nodal  points  of  these  elemental  areas.  The 
numerical  procedure  is  based  on  the  use  of  Eq.s(C7-94)  with  y^  =  0  for 

purposes  of  determination  of  and  its  derivatives  at  points  on  the  free 

*  ^ 

surface.  The  function  f  is  then  evaluated  at  the  nodal  points  by  use  of 


Eq. (50) . 

The  solution  for  the  second-order  source  strength  is  obtained  numerically 
through  the  integral  equation,  E q . ( 79 ) .  That  is,  Eq . ( 79 )  is  written  as 

J  »/,  V  * 

.  K 


£.  =  hLi  torik(*v)  f’ 


(95/ 


where  &LJ  (+v)  is  defined  by  replacing  V  by  in  Eq .  ( S3 )  and 


hit  -  t  i nti  ccsh^d(h*yt)]]3 


zaxk 


(96  1 


Once  the  CK  -matrices  in  (05)  have  been  evaluated  it  may  be  solved  for  the 
source  strength  distribution  on  the  immersed  surface. 


3  1 


* 


The  hydrodynamic  pressure  and  resultinq  forces  carried  to  the  second 
order  as  given  in  Eq.(59)  and  (69),  respectively,  require  only  the  second- 
order  potential  rather  than  its  derivatives.  Once  f^  is  known,  u|  roay  be 
obtained  through  Eq.(73)  which,  when  written  in  indicial  form,  becomes 


~  A* ( *v>  6* 


(97) 


where  fiuC+v)  i s  defined  by  Eq. (91 ) . 


NUMERICAL  RESULTS 

A  computer  code  has  been  developed  for  the  special  case  of  a  vertical 

circular  cylinder  (pile).  MacCamy  and  Fuch's  solution  was  used  for  the 

* 

first-order  solution  and  for  purposes  of  evaluating  f  .  The  second-order 
solution  was  calculated  on  the  basis  of  the  distributed  singularities  as 
qi ven  by  Eq. ( 73) . 

The  particular  example  calculation  was  carried  out  for  a  cylinder  nlaced 
in  water  one  radii  deep.  The  dimensionless  frequency  parameter  was  varied 
over  the  range  C<.  V <  1.2.  The  results  for  the  first-order  horizontal  force 
coefficient,  the  second-order  horizontal  force  coefficient,  steady-state 
force  coefficient  and  phase  anqles  is  shown  in  Fiaure  3.  The  results  appear 
to  indicate  that  the  second-order  effects  are  greatest  at  the  lower  freq¬ 
uency  ranqe.  i.e.,  for 

In  Figure  A  and  5  the  dimensionless  horizontal  force  is  plotted  for  a 
complete  wave  cycle  for  two  different  frequencies  (wave  lenqths).  These 
results  which  correspond  to  rather  steep  wave  show  a  rather  sizeable  non¬ 
linear  effect.  Also,  the  effect  of  including  second-order  effects  is  to 
shift  the  maximum  force  towards  the  phase  of  the  wave  crest. 
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